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The employment of path entangled multiphoton states enables measurement of phase with en- 
hanced precision. It is common practice to demonstrate the unique properties of such quantum 
states by measuring super-resolving oscillations in the coincidence rate of a Mach-Zehnder interfer- 
ometer. Similar oscillations, however, have also been demonstrated in various configurations using 
classical light only; making it unclear what, if any, are the classical limits of this phenomenon. 
Here we derive a classical bound for the visibility of super-resolving oscillations in a Mach-Zehnder. 
This provides an easy to apply, fundamental test of non-classicality. We apply this test to experi- 
mental multiphoton coincidence measurements obtained using photon number resolving detectors. 
Mach-Zehnder super-resolution is found to be a highly distinctive quantum effect. 

PACS numbers: 42.50.-p,42.50.Xa,03.65.Ud 



Inroduction. — Experimental observations which ex- 
hibit uniquely quantum mechanical behavior lie at the 
heart of quantum optics. Such observations establish 
the way we difFrentiate between the quantum and classi- 
cal worlds. Notable examples include Bell's inequalities 
[H , photon anti-bunching [3] and sub-Poissonian number 
statistics 0. Negativity of the Wigner function is also 
widely used to verify non-classicality ^4j] . The derivation 
of additional tests involving non-classical behavior has 
motivated many studies |5l4lG|. 

Recently, there has been a surge of interest in path 
entangled two- mode states of the form |iV :: 0)a b = 
^ (1^, 0)a.b + |0, N)a.,b) known as NOON states [Hi. In 
such states, all N photons acquire phase as a collec- 
tive entity thus enabling phase super-sensitivity and 
sub- wavelength lithography A number of experi- 

mental implementations of NOON states exist to date 
In these experiments the high NOON state 
fidelity was demonstrated by measuring a sinusoidally 
varying N photon coincidence rate with a reduced pe- 
riod corresponding to a wavelength X/N. Such behavior, 
coined 'phase super-resolution' [ll] is generally accepted 
as a signature of NOON states. However, it has been 
shown that similar measurements exhibiting high visibil- 
ity super-resolution can be obtained using classical light 
only. This has been demonstrated using multiport lin- 
ear interferometers with coincidence counting il7J o r in a 
scheme for nonlinear sub-Rayleigh lithography [iSl [l9| . 

It is therefore natural to ask whether super-resolution 
can be used at all as a direct, non- ambiguous test for 
non-classicality. Here we answer this question in the af- 
firmative. We consider a Mach-Zehnder geometry simi- 
lar to the one employed in the initial proposal for quan- 
tum lithography [12| and subsequently used in all NOON 
state generation experiments (see Fig. [IJ . We show that 
for any classical state with a positive-definite well be- 
haved Glauber-Sudarshan P representation [2^, there is 
a strict upper limit for the visibility of super-resolving 
fringes. The bound is derived analytically for arbitrary 



coincidences with m, n photons in the two output ports 
respectively. We apply the bound to experimental mul- 
tiphoton coincidences obtained using both classical and 
quantum light. The classical results are found to ap- 
proach the bound from below. For the quantum case, 
we employ NOON states obtained in an experiment we 
have recently reported 15| and a violation of the classical 
bound with states of up to five photons is demonstrated. 

Derivation. — Consider a balanced Mach-Zehnder In- 
terferometer (MZI) whose input modes are denoted a, b 
(see Fig. [T]). The coincidence rate with m,n photons in 
detectors Di,D2 respectively is denoted Cm,n{<f) where 
the total number of photons is TV = n+m. When employ- 
ing a NOON state input, Cm.n exhibits perfect sinusoidal 
fringes of the form Cm.nif) c>c 1 -I- cos{N(p ~ S) where 
i5 = 0, TT depending on the values of m, n. This pattern 
exhibits A^-fold phase super-resolution with 100% visibil- 
ity. Our goal is to derive a simple unambiguous criterion 
based on the visibility of Cm.n for distinguishing between 
quantum and classical input states. To this end, consider 
the Fourier series of an arbitrary coincidence function. 



k=0 



Ak cos(fc(p - Sk). 



(1) 



We are interested in the visibility of the A^-fold oscilla- 
tions, |AAr/j4o|, which will be referred to as the TV-fold 
visibility. For pure sinusoidal oscillations with a constant 
background the iV-fold visibility coincides with the con- 
ventional definition of visibility. Considering the A'^-fold 
visibility allows application of this bound to situations in 
which the measured coincidence has contributions oscil- 



lating at a number of different frequencies (see [15| for 
example). The m,n coincidence rate for a MZI is given 
by [3, 



C'm.n {V) = Tr 



(2) 



where U {f) is a unitary operator describing the MZI 
p3j . pa,b = \i^)a,b{''P\ is the input state density ma- 
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FIG. 1: (color online) Experimental setup, (a) Schematic of 
the setup, a Mach-Zehnder interferometer with photon num- 
ber resolving detectors at the output ports, (b) The actual 
setup used for obtaining the classical light results shown in 
Fig. [3] A horizontally (H) polarized pulsed Ti:Sapphire os- 
cillator with 120 fs pulses @1 MHz repetition rate (reduced 
from 80MHz using a pulse picker) is used as the light source. 
The MZI is implemented in a coUinear inherently phase stable 
design using polarization. The phase, ^ is applied using a liq- 
uid crystal phase modulator at 45 deg. (X, Y polarizations). 
The photon number resolving detectors are Hamamatsu multi 
pixel photon counters (MPPC) module C10507-11-050U 
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FIG. 2: The classical visibility bound, F^'" (Eq. (gj) as a 
function oi N = m + n. For each A^, the bound is shown for 
two choices of m, n: (i) m,n — [A/2J , [A/2] i.e. the detection 
events are divided as equally as possible between Di and D2 
(see Fig. [l]), (ii) m,n — N,0 i.e. all N photons are detected in 
Di and zero are detected in D2 . For a given N, the classically 
obtainable visibility for case (i) is always higher than for case 
(ii). However, for both cases F^'" decays exponentially with 
N. 



trix and tt^*-^^ are the positiv e op erator valued measures 
(POVM's) of detectors L'i(2) uM- For the present deriva- 
tion we consider ideal photon number resolving detectors 
implying ■ir^ = tt^ = \n){n\. First, we derive the A^-fold 
visibility for an input state of the form \ip)a,b = \ct)\0)a,b- 
Here, \a) is an arbitrary non- vacuum coherent state de- 
fined \a) = Er=oe"^'"''7=Tl")- Substitution in Eq. (g]) 
yields. 



1 ' m ' 



|a|2(™+")| sin((p/2)|2'"| cos(¥>/2)|^ 



2{m+n) 



(3) 



We denote the A^-fold visibility for each photon number 
pair, m, n by F™'". Using simple combinatorial consid- 
erations we obtain 



r=0 ^ ^ 



2™ 

n + m — r 



(4) 



We note that F^^^'" 



does not depend on the specific choice 
of a. In the following we prove that F™'" is, in fact, the 
largest A^-fold visibility obtainable for an arbitrary classi- 
cal input state. We now consider an input state consist- 
ing of two arbitrary coherent states |i/')a,b = |a)|/?)a,fc- 
We denote the A^th fourier component of Cm,n{f) for 
this state by Using some straightforward but 

somewhat lengthy algebra it can be shown that for such 
an input state the TV-fold visibility is always bounded by 



I.e. 



^Q./3,m.n 



IK 



(5) 



Finally, we prove the bound holds for general clas- 
sical inputs of the form p^'^ = / P(a, /3)a_f,|a)a(a| <8> 
/3)f,(/3|d^a(i^/3 with a positive-definite, well behaved P 
function P(a,/3)a,b. These are precisely the quantum 
mechanical states which posses classical analogues [2^ . 
Formally the iV-fold visibility is a functional of P which 
we denote Vm,n\P{oL,^)\. The bound ensues from the 
following inequalities. 



T4.,„[P(a,/3)] = 



< 



< 



/P(a,/3)Ao" 



/P(a,/3) 



a,/3,m.n 



(6) 



A, 







JP{a,(3)A'^'^'"''''d^ad^(3 



/P(a, /3)F™'M; 



m,n ^a,/3,m,n 72 







d'ad^(3 



JP{a,(3)A^'^'"'^''d^ad^(3 



= F" 



using the positive definiteness of P and Eq. 



Note 



that Aq''^'"^'"' is positive by definition. This completes 
the proof and the main result of this letter. 

The bound, F^'", is plotted in Fig. [2]for a selection of 
TO, n pairs. In addition, the numerical values for N < 5 
are given in Table H] The function F^'" decays expo- 
nentially with N = m + n. However, for a given N 
the obtainable visibility is much higher when the pho- 
tons are distributed equally (or almost equally for odd 
A^) between the two output modes. Interestingly, our 
result for pairs of the form TV, coincide with those ob- 
tained by Bentley and Boyd in a proposal for sub- 
wavelength lithography using classical light thus showing 
their scheme reaches the classical bound. 

We note that although the derivation uses a single 
spatial and spectral mode, its generalization to multi- 
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TABLE I: Classical visibility bound, 
cal values for N < 5. 



Eq. (g]). Numeri- 



K 2k 
MZ phase [rad] 



ji 2 II 
MZ phase [rad] 



FIG. 3: (color online) Experimental coincidence scans using 
classical light, (a-d) Coincidence rates (raw data, normalized) 
as a function of Mach-Zehnder phase with m 'clicks' in Di 
and n 'clicks' in D2 for various values of m,n (see Fig[l]). 
The total number of photons m + n is denoted by A*'. The 
input state is \i')a,b ~ \ct}\0)a,b where |a) is a coherent state. 
The solid lines are obtained using an analytical model of the 
detector POVM's accounting for cross-talk, dark counts and 
losses [23] . Error bars are smaller than the displayed markers. 
(a.) N ^ 2 (i) m,n = 2,0, (ii) m,n = 1,1. (b) = 3 (i) 
m,n = 3,0, (ii) m,n = 2,1. (c) A^ = 4 (i) m,n = 4,0, 
(ii) m,n — 2,2, (iii) m,n = 3,1. (d) N=5 (i) m,n — 5,0, 
(ii) m,n — 3,2, (iii) m,n — 4,1. The classical visibilities in 
Fig. |4)ii (blue bars) are obtained using these measurements, 
(e) Three-fold sinusoidal oscillations used for obtaining the 
experimental visibility of the m, n = 2, 1 classical case. To 
obtain the black squares, three copies of the raw data shown in 
panel bii were shifted by 0, ^ and superimposed. The red 
line is the fitted Fourier series. The obtained 3-fold visibility is 
45.22%. (f) Four-fold sinusoidal oscillations used for obtaining 
the experimental visibility of the m,n — 2,2 classical case. 
Four copies of the m,n = 2,2 data shown in panel cii were 
shifted by 0, |-,7r, ^ and superimposed. The obtained 4-fold 
visibility is 29.19%. 



pie modes is straightforward [25j and does not alter the 
final result. We have implicitly assumed the use of broad- 
band detectors such as those typically employed for single 
photon detection. If one allows use of narrow-band mul- 
tiphoton detection (such as atomic two photon absorp- 
tion) and appropriate engineering of the input state then 
higher visibilities may be obtained [Hill. The classical 
visibility limits of multiphoton interferences have been 
previously discussed 



28 



|29|. However, these analyses 
were limited to input states with randomly fluctuating 
phases and considered the standard rather than iV-fold 
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visibility which in our case can always be 100% (as seen 
experimentally in Fig. [3|). We point out that phase 
sensitivity beyond the standard quantum limit (known 
as phase super-sensitivity) is a distinct quantum effect 
and can be exploited for quantum noise reduction. How- 
ever, it cannot be inferred directly from the visibility 
and requires additional knowledge about the efficiency 



171 121| | , making it much less straightforward as a test for 
non-classicality. The sensitivity criterion is also typically 
harder to beat experimentally than the present bound. 

Experimental results. — We now demonstrate the appli- 
cability of the bound using experiments done in both the 
classical and quantum regimes. In the classical case, the 
bound is approached from below by using a MZI with a 
coherent state input of the form |'(/')a,b = |Q;)|0)a,6. The 
MZI is implemented in a polarization based geometry as 
shown in Fig. [TJj. We measured experimental multipho- 
ton coincidences using two multi-pixel photon number re- 
solving detectors . The coincidence sweeps are shown 
in Fig. [3^-d. For a given N it is readily observable that 
the features are widest when the photons are all detected 
in one output. It is therefore not surprising that the vis- 
ibility bound is lowest in this case. The A^-fold visibil- 
ity was obtained by fitting the experimental curves to a 
Fourier series of the form in Eq. ([T]) truncated at fc = iV. 
For a given N, the fit was applied to N copies of the data 
shifted by 0, 27r-^, . . . 27r-^^^ and superimposed (see Fig. 
[5j3-f). This eliminates the lower frequencies without af- 
fecting the A'^-fold visibility thus improving the accuracy 
of the numerical fit. Although in this case the slowly os- 
cillating terms were eliminated in post-processing, it can 
be shown that for any value of m, n there exists a classi- 
cal input state which yields a coincidence rate of the form 
Cm,n{'p) oc 1 -I- r™'" cos{Nip + S) i.e. pure A^-fold oscil- 
lations with a constant background. The experimentally 
obtained visibilities are shown in Fig. |4jii (blue bars). 
The measured visibility is lower than the classical bound 
for all TO, n as predicted. We note that limited overall ef- 
ficiency doesn't affect the classically obtainable visibility 
since coherent states are transformed to coherent states 
by loss '30]. Additional experimental imperfections such 
as cross-talk and dark-counts only serve to reduce the 
A^-fold visibility. Thus, the bound is expected to hold for 
raw experimental data and doesn't require compensation 
for detector imperfections. 

Surpassing the bound requires use of non-classical 
light. To demonstrate this we use path-entangled mul- 
tiphoton states with high NOON state fidelity that we 
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FIG. 4: (color online) Comparison between the classical 
bound and experimental visibilities. Visibilities as a function 
of the number of coincidence 'clicks' in each of the MZI out- 
puts, (i) The classical bound, F^'", Eq. Q. (ii) Experimen- 
tal visibilities obtained with NOON states using non-classical 
light [15]. All the visibilities violate the classical bound (ex- 
cept 4,0 which was not measured and 1, 1 which cannot be 
surpassed). Error bars indicate ±o" statistical uncertainty, 
(iii) Experimental visibilities obtained using classical inputs 
(error bars are negligible). The experimental visibilities were 
obtained from the coincidence measurements in Fig. O These 
visibilities are lower than the classical bound as expected, the 
discrepancy being mainly due to the cross-talk between de- 
tector pixels. 



have recently reported [l5| . Briefly, the experimental ge- 
ometry was similar to the one used in Fig. [TJ however 
one of the inputs was non-classical resulting in a highly 
entangled state in modes c, d of the MZI. The visibilities, 
obtained from raw data, are in violation of the classical 
bound. This verifies the quantum nature of the results 
and demonstrates the feasibility of surpassing the clas- 
sical bound for a wide range of photon numbers in a 
realistic experimental setting, see Fig. |4ji (red bars). 

Conclusion. — We have derived a simple, practical test 
for non-classicality using MZI coincidence measurements 
and demonstrated its use. The quantity of interest is the 
visibility of A^-fold oscillations in the N photon coinci- 
dence rate. The classical bound on this quantity drops 
exponentially to zero with the total photon number N 
while the quantum limit remains constant. This makes it 
an excellent indicator of non-classicality. It can be shown 
that only the NOON component of a given state in the 
MZI contributes to the A^-fold oscillatory term. Thus, 
our results reflect the fact that the NOON state overlap 
resulting from any classical input state becomes exceed- 
ingly low as N grows. We expect this result to serve as 
a benchmark in future experimental demonstrations of 



super-resolution. 
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